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Abstract 



We study a random positive definite symmetric matrix distributed according to a real 
Wishart distribution. We compute general moments of the random matrix and of its inverse 
explicitly. To do so, we employ the orthogonal Weingarten function, which was recently 
introduced in the study for Haar-distributed orthogonal matrices. As applications, we give 
formulas for moments of traces of a Wishart matrix and its inverse. 



1 Introduction 

1.1 Wishart distributions 

Let d be a positive integer. Let Sym((i) be the M-linear space of d x d real symmetric matrices, 
and $7 = Sym^{d) the open convex cone of all positive definite matrices in Sym((i). Let a = 

{(^ij)i<i,j<d G ^, and let 

^ [12 d-l] fd-1 

Then there exists a probability measure Wd,i3,a on such that its moment generating function 
(or its Laplace transform) is given by 



Jn 



where 6 is any d x d symmetric matrix such that a^^ —6^0,. We call Wd,/3.a the real Wishart 
distribution on il. with parameters a) . 

We call a random matrix G a real Wishart matrix associated with parameters (/3, a) 
and write W ~ a; M) if its distribution is 2Bd,/3,o-. Thus the moment generating function 

for W is given by 

E[e*'^(^^)] =det(/d-^a)-^ 
with 9 being as above. Here E stands for the average. 



^ MSG 2000 subject classifications: Primary 15A52; secondary 60E05, 20C30. 
Keywords: Wishart distribution, Gelfand pair, zonal polynomial, orthogonal group, hafnian. 
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If 2/3 is a positive integer, p = 2(3 say, then a Wishart matrix is expressed as follows. Let 
Xi, . . . , Xp be d-dimensional random column vectors distributed independently according to the 
Gaussian distribution Nd{'^,^a). Then 

W = XiX\ + ■■■ + XpX'p, 

where X^ is the transpose of Xi, i.e., a row vector. 

If /3 > (not necessarily an integer), the distribution Wd,i3,a has the expression 

WdM^w) = f{w;d,(3,a)£{dw), 

where f{w;d, (3,a) is the density function given by 

(1.1) f{w;d,^,a) = rrf(/3)-i(deta)-^(detw;)''-'^e-*'^(""'^) {w G O) 

with the multivariate gamma function 

r,(/3) = 7r'^('^-^)/^nr(/3-i(j-i)). 

Here £ is the Lebesgue measure on Sym((i) defined by 

2{dw) = dujjj with w = {wij)i<ij<d. 

^<i<j<d 

Likewise, a complex Wishart distribution is defined on the set of all d x d positive definite 
hermitian complex matrices. Given a Wishart matrix W, the distribution of the inverse matrix 
W^^ is called the inverse (or inverted) Wishart distribution. We denote by Wij and W^^ the 
(i,j)-entry of W and W''^, respectively. 

The Wishart distributions are fundamental distributions in multivariate statistical analysis. 
We refer to jMuj . The structure of Wishart distributions has been studied for a long time, 
nevertheless, a lot of results are obtained only recently. We are interested in moments of the 
forms E[P(1^)] and E[P(T4^~-'^)], where P{A) is a polynomial in entries Aij of a matrix A. In 
particular, we would like to compute general moments 

nWi.j, Wi,j, ■ ■ ■ Wi,j,] and E[W''^^ W''^' ■ ■ ■ W'^^^] 

for W and W~^, respectively. 

Von Rosen |Voj computed the general moments of low orders for W^^. Lu and Richards 
|LR] gave formulas for W by applying MacMahon's master theorem. Graczyk et al. [GLMl] 
gave formulas for W^^ in the complex case by using representation theory of symmetric groups, 
while they [GLM^ gave results for only W (not VF~^) in the real case by using representation 
theory of hyperoctahedral groups. Letac and Massam |LMlj computed the moments E[P(VF)] 
and E[P(IF~^)] in both real and complex cases, where the P is a polynomial depending only on 
the eigenvalues of a matrix. Furthermore, a noncentral Wishart distribution was also studied, 
see [LM2] and [KNT]. 
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1.2 Results 

Our main purpose in the present paper is to compute the general moment 

E[W'^^^W'^^^ ■ ■ ■ W''''^'^] 

for an inverse real Wishart matrix = (W^^). As we described, in the complex case Graczyk 
et al. [GLMlj obtained formulas for such a moment by the represention-theoretic approach. 
Our main results are precisely their counterparts for the real case, which had been unsolved. 

To describe our main result, we recall perfect matchings. Let n be a positive integer and put 
[n] = {1, 2, . . . , n}. A perfect matching m on the 2n-set [2n] is an unordered pairing of letters 
1,2,..., 2n. Denote by A4{2n) the set of all such perfect matchings. For example, A4(4) consists 
of three elements: {{1, 2}, {3, 4}}, {{1, 3}, {2, 4}}, and {{1, 4}, {2, 3}}. 

Given a perfect matching m G M.{2n), we attach a (undirected) graph G = G{m) defined as 
follows. The vertex set of G is [2n]. The edge set of G is 

{{2k -1,2k} I ke[n]}u{{p,q} \ {p,q}£m}. 

Then each vertex has just two edges, and each connected component of G has even vertices. 
We denote by fi;(m) the number of connected components in G(m). For example, given m = 
{{1, 3}, {2, 7}, {4, 8}, {5, 6}} G A^(8), the graph G(m) has two connected components (where 
one component has vertices 1, 2, 3, 4, 7, 8 and the other has 5, 6), and therefore K(m) = 2. 

Now we give a formula of the general moments for W. It would be reader-friendly to describe 
the results of both the real Wishart law itself (Theorem[T]) and its inverse (Theorem[2|) on display. 

Theorem 1. Let W = (W^jj)i<jj<d ~ Wd{(3, a; M). Given indices ki,k2, ■ ■ ■ , /c2n. from {1, . . . , d}, 
we have 

(1.2) ElWk,k,Wk,k,---Wk,„_,k,J=2-^ (2/3)"^"^ n ^fep*^.- 

meM{2n) {p,q}em 

For example, since «;({{1, 2}, {3, 4}}) = 2 and 3}, {2, 4}}) = 4}, {2, 3}}) = 1 we 

have 

Theorem [T] is not new. Indeed, it is equivalent to Theorem 10 in |GLM2j . Moreover, Kuriki 
and Numata [KNlj extended it to noncentral Wishart distributions very recently. We revisit it 
in the framework of alpha-hafrians. We develop a theory of the alpha-hafnians in Section [2l and 
apply it to the proof of Theorem [1] in Section [3l 

The following is our main result. Let o"*-^ be the (i, j)-entry of the inverse matrix cr~^. 

Theorem 2. Let W ~ VFrf(/5, cr; M). Put ^ = j3 — and suppose 7 > n — 1. Given indices 
ki,k2, ■ ■ ■ , k2n from {1, . . . , d}, we have 

(1.4) ]E[vi/'^'ifc2|^fe3fc4 ...p^fc2n-ifc2n] ^ ^ Wg(m; 7) 0"*=" . 

meA4(2n) {p,9}Gm 

Here Wg(m;7) is defined in Section\^ below. 
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For example, for 7 > 1 we will see 

27 - 1 



Wg({{l,2},{3,4}};7) 



"7(7-l)(27 + l)' 
Wg({{l, 3}, {2, 4}}; 7) = Wg({{l, 4}, {2, 3}}; 7) =-^---l__^, 

and therefore we have 

jg|-^fcifc2^A:3fc4j _ \ (^^2j — I^(jkik2^k3k4 _|_ ^kik-i ^k2k4, _|_ ^kik4^k2k3y 

The quantity Wg(m; 7) is a slight deformation of the orthogonal Weingarten function. The 
function was introduced by Collins and his coauthors \CM\ ICS] , in order to compute the general 
moments for a Haar-distributed orthogonal matrix. In general, Wg(m;7) (m G A4{2n)) is 
given by a sum over partitions of n, and derived from the harmonic analysis of the Gelfand pair 
{S2n,Hn), where S2n is the symmetric group and Hn is the hyperoctahedral group. Amazingly, 
the same function thus appears in two different random matrix systems. In Section [H we review 
the theory of the Weingarten function developed in \CM\ IMat2j . and, in Section [5l we prove 
Theorem [2j 

In Section [6l we give applications of Theorem [1] and Theorem [2j In particular, we obtain 
results of Letac and Massam ^jMJJ as corollaries of Theorem [1] and Theorem [21 

2 Alpha-hafnians 

2.1 An expansion formula for alpha-hafnians 

Let A be a 2n X 2n symmetric matrix A = (^pq)p,<je[2n]- Let a be a complex number. We define 
an a-hafnian of A (see [KN2j ) by 

hf„(^)= Yi n ^pi- 

meM{2n) {p,q}em 

The ordinary hafnian of A is nothing but hfi(74). For example, if n = 2, 

Ma{A) = a^AuAs4 + 0^13^24 + 0^14^23- 

We remark that hfQ,(^) does not depend on diagonal entries Aii,A22, • • • , ^2n,2n- Note that the 
right hand side in (fL2]) is equal to 2~"hf2/3(crfcpfcjp gg[2„]. 

Proposition 1. Let A = (^pg)p,gg[2n] be a symmetric matrix. Let D = (^pg)p,ge[2n-2] ■ Por each 
j = 1,2, . . . ,2n — 2, let B^^^ he the symmetric matrix obtained by replacing the jth row/column 
of D by the (2n — l)th row/column of A. In formulas, B^^^ = {Bpq)p q^y2n-2\ given by 



PI 



A2n-i,2n-i if P = j and q = j, 

A2n-i,q ifp = j andq^j, 

Ap^2n-i ifPT^j and q = j, 

^P,q ifPT^j andq^j. 
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Then we have 

2n-2 

(2.1) hf„(A) = ^ 2nM„(5(-')) + a^2n-l,2nM,(I?). 

We call ()2.ip an expansion formula for an a-hafnian with respect to the (2n)th row/column. 

Proof. For each j = 1, 2, . . . , 2n — 1, we set 

Mj{2n) = {me M{2n) | {j, 2n} G m}. 

Then M{2n) = lJj=T^ A^j(2n). We define a one-to-one map m n from Mj{2n) to A^(2n - 2) 
as follows. 

First, suppose j = 2n — 1. Given m G M2n-i{2n), we let n G M{2n — 2) to be the perfect 
matching obtained from m by removing the block {2n — l,2n}. It is clear that the mapping 
M2n-i{2n) 9 m n G M{2n — 2) is bijective and that K(m) = K(n) + 1. 

Next, suppose j G [2n — 2]. Given m G Mj{2n), we let n G M{2n — 2) to be obtained by 
removing the block {j, 2n} and a block {i, 2n — 1} (with some i G [2n — 2]) and by adding {i, j}. 
It is easy to see that this mapping M.j{2n) 9 m i-)- n G M{2n — 2) is bijective, K(m) = ^(n), and 

For example, consider m = {{1, 4}, {2, 5}, {3, 6}}. Then m G A^3(6), and we obtain n = 
{{1,4}, {2, 3}} G M(4). Therefore, we have ^(m) = 1 = ^(n) and Y\{p,q}&m^pq = ^14^25^36 = 

A r(3) r(3) _ 4 n r(3) 

Using the correspondence M.j{2n) 9 m o n G M.{2n — 2) with j = 1, 2, . . . , 2n — 1, it follows 
that 

2n-l 

hf«(^) = ^ A,- 2n a^^^^ J] 

i=l meA1j(2n) {p,g}Gm 

{p,g}7^{j,2n} 

=^2n-l,2n ^ a''(")+^ J] ^f'' 
nGA4{2n-2) {p,g}en 

2n-2 

J=l neA1(2n-2) {p.gjen 

which is equal to A2n-i,2nO: ■ MaiD) + X:?=T' ^j,2nhf„(i?(^)). □ 

2.2 Another expression for «-hafnians 

Let Sn be the symmetric group on [n]. Each permutation tt is uniquely decomposed into a 
product of cycles. For example, '?i" = (56i432)^'^6is expressed as vr = (1 ^ 5 — )• 3 ^ 1)(2 — > 
6 —7- 2) (4 ^ 4). Denote by C{tt) the set of all cycles of vr, and let i^iir) be the number of cycles 
of vr: i/(7r) = |C7(7r)|. 
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Let A = {Apq)pq^^2n] be a symmetric matrix. For each k,l G [n], we denote by ^[A;, Z] the 
2x2 matrix 

A[k,l] = (^f-^'^'-' ^f-^'A. 

V A2k,2l-1 A2k,2l J 

For a cycle c = (c^ —>■ ci —>■ C2 —>■•••—>■ c^) on {1, ... , n}, we put 

Pc{A) = tr {A[ci,C2]JA[c2,C3]J ■ ■ ■ A[cr, ci] J), with J = 

In particular, P(^ci^ci){^) = tr (^[ci, ci] J) = 2A2ci-i,2ci for a 1-cycle (ci — > ci). It is easy to 
see that Pc{A) can be written 

summed over (j2fe-i)j2fc) £ {(2cfc — l,2cfc), (2cfe,2cfc — 1)} (fc = 1,2, . . . ,r). For a permutation 
TT G -Sn, we define 

ceC(7r) 

Similarly, given an r-cycle c = (c^ ^ ci ^ C2 c^), we let c,. to be the largest number 

among {ci, C2, . . . , Cr}. We define Qc{A) as follows: If r = 1 then Qc{A) = A2ci-i,2ci', if r > 2 
then 

Qc{A) = ^ • • • ^ A2c^-lJiAj^j^Aj^j^ ■ ■ ■ Aj2^_2,2Cr.) 

(ilj2) (i2r-3,i2r-2) 

summed over {j2k-i,j2k) G {(2cfc - 1, 2cfc), (2cfc, 2cfc - 1)} (fc = 1, 2, . . . , r - 1). As P7r(-4), we 
define 

QAA)= H Q,{A). 

ceC(7r) 

For example, for a cycle (3 ^ 2 — )■ 1 ^ 3), we have 

Qc{A) = ^ ^ A^jj^Aj^j^Aj^Q 

(iiJ2)e{(3,4),(4,3)}03j4)e{(l,2),(2,l)} 
=^53^41^26 + ^54^31-426 + ^53^42^16 + ^54^32^16- 

Lemma 2. Lei c = (c^ ci —)■ C2 c^) 6e a cycle. Then 

PM) = Qc{A) + Q,-i{A), 
where = {cr ^ ■ ■ ■ ^ C2 ^ ci ^ Cr). 

Proof. Suppose Cr is the largest number in {ci, . . . , c^}. We can express 

Pc{A) = ^2 ^2Cr-lJl^J2j3 ■ ■ ■ ^j2r-2,2CT. + ^ ^2Cr ,J1 ^j2,J3 ' ' ' ^j2r-2,2c,.-l ) 

,32 , ■ ■ ■ , J2r- - 2 il ,i2 , ■ ■ ■ :i2r - 2 

summed over (j2fe-i, i2fe) G {(2cfc - 1, 2ck), {2ck,2ck - 1)} (fc = 1, 2, . . . , r - 1). Here the first 
sum coincides with Qc{A), while the second one does with Q^-i{A). □ 
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Proposition 3. Let A = (^pij)p,q6[2n] ^6 o, symmetric matrix. Then 



This is a key lemma in our proofs of Theorem [T] and Theorem [2j We show this proposition 
in the next subsection. 

Remark 1. Let A = (Ajj)i<j j<„ be a complex matrix and a a complex number. An a- 
permanent of A is defined by 

n 

perJA)= a^^-^n^-W- 
It intertwines the permanent and determinant: 

n 

peri(^) = per(yl) = ^ JJ^^-W ^nd per_i(^) = (-l)"det(A). 



It is also called an a -determinant. See |Vej and also [Sh] , Alpha-hafnians are generalizations of 
the alpha-permanents in the following sense. Given a matrix A = (vljj)i<i,j<n, we define the 
2n X 2n symmetric matrix B = (i3pg)i<p.g<2n by 

B2i-i,2j~i = B2i,2j = and B2i^i,2j = ^2i,2i-i = Aij for all i, j = 1, 2, . . . , n. 

Then, since Qc{B) — ^(;r,ci^ci,c2 • • • ^cr-i,cr foi' c — (Cf. — ^ ci — y C2 — y ■ ■ ■ — y Cf ), it follows from 
Proposition [3] that hia{B) = perQ,(A). Thus any a-permanent can be given by an a-hafnian. 

Remark 2. Let B = (-Bpij)p,ge[2n] be a skew-symmetric matrix and let a be a complex number. 
In [Matlj . an a-pfaffian of B was defined. In a similar way to the proof of Proposition [31 we can 
see that the definition in [Matl] is equivalent to the expression 

pUB)= Y1 (-«)"^"^sgn(m) n 

m&M{2n) {p,q}€M{2n) 

Here, for m = {{m(l), m(2)}, . . . , {m(2n — l),m(2n)}} we define 

sgn(m) n ^^"^ = ^8^(^(1) m(2) m(2n)) " " " " ^' 

{p,q}&M{2n) ^ ^ ' ^ 



m(l)m(2) ■ ■ ■ -Dm(2?i-l)m(2n)- 



When a = — 1, the a-pfaffian is exactly the ordinary pfaffian. Moreover, as a-hafnians are so, 
the a-pfaffians are generalizations of a-permanents. 
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2.3 Proof of Proposition [3] 

Put 

TT^Sn Tr(zSn 

for any n > 1 and any symmetric matrix A of size 2n. Here the second equality follows from 
Lemma [2j 

Let 5W,5(2),...,5(2n-2)^i^ be as in Proposition [TJ In order to obtain Proposition [21 it is 
enough to show the recurrence formula 

2n-2 

(2.3) M^iA) = A,-2„hf,(5(^)) + aA2n-l,2n^UD). 

To see ()2.3p . we will show a recurrence formula involving Qc{A) and Pc{A). For each A; G [n], 
we denote by S^'^ the subset of permutations in Sn such that 7r{k) = n. Note Sn = |Jfc=i 

Let A; € ['^ — 1] and let vr € S'l'^^ Let Un(7r) S ^^(Tr) be the cycle including the letter n, which 
is of the form 

UniT^) = (n — Ci — C2 — 5- • • • — > — )■ — 5- n) , 
with (possibly empty) distinct ci, . . . , Cj. € [n] \ {/c, n}. Then, define 

UniT^) = {n ^ Cr ^ ■ ■ ■ ^ C2 ^ Ci ^ k ^ n), 

and let tt be the permutation obtained by replacing Uniir) in vr by Un{iT). Note that nn(7r) = 
u„(7r), and that tt = vr if and only if Un{T^) is a 2 or 3-cycle. The map vr i— )• vr is an involution 
on Sn^. For example, given vr = (7 3 ^ 1 2 ^> 7)(6 4 ^ 6)(5 ^ 5) e 5/, we have 
# = (7^1^3^2^ 7)(6 ^ 4 ^ 6)(5 ^ 5). 

In general, for the cycle n„(vr) = (n — )■ ci — )■ C2 — > • • • ^ — )• A: — )• n) with k ^ n, we see 
that 

Qu,i{tt) -l,Ji^J2,j3 ' ' ' ^j2r,2k-lA2k,2n + ^2n-l,ji^j2,i3 ' ' ' ^j2r,2A:^2fc-l,2n) 

(jlj2) {j2r-lj2r) 

_ \^ fry{2k) r.{2k) „(2fc) . „{2fc-l) „{2fc-l) R(2fc-1) . N 

(il,i2) {J2r-l,j2r) 

summed over 

{hp-iihp) e {(2cp - l,2cp), (2cp,2cp - 1)} (p = 1,2, . . . ,r). 

Similarly, 

Qu„(7r)(^) = E "' E (^2n-lj2r ■ ■ ■ ^i3,i2^jl,2fc-1^2fc,2n + ^2n-l,j2r ■ ■ ■ ^i3,i2^Ji,2fc^2fe-l,2n) 

(ilJ2) {i2r-lj2r) 

_ \^ CR(2fc) „(2fe) „(2fc) . „(2fc-l) „(2fc-l) „(2fc-l) , ^ 

- 2^ ■■■ 2^ ^-^2fcj2r ■ ■ ■ ^i3 J2-^jl,2fe-1^2fc,2n + i^sfc-ljar ■ ■ ■ -'^j3J2 ^h,2k ^2k-l,2n) 

{31,32) {i2r-l,i2r) 
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Therefore, we have 

(2.4) Qn„w(^) + Qu^nM) = A2fc,2n^<w(i?('')) + A2fc_i,2„P< (.) (i^^''-^^). 

Here tiJi(vr) is the cycle obtained from tin(vr) by removing the letter n: u^(vr) = (A; — )• ci ^ C2 — > 
■ ■ ■ — > Cr ^ /c). We note that the mapping vr vr' := ^n(7r) nc6(7(7r)\{M„(7r)} bijective 

map from Sn^ to Sn-i, and that z^(7r) = i^{tt'). 

Now we go back to the proof of (|2.3p . We rewrite 

n-l 

The first sum is equal to 

a'^(-')+'g.'(^)Q(n)(^) = aA2„_l,2nhfa(^) 

by a natural bijective map 5^"^ — )• while, since the map vr i-^ vr is bijective on each s!^^ , 

the terms corresponding to /c € [n — 1] in the second sum are equal to 

vresi") cGC(7r)\{«„(7r)} 



(fc) ceC{7r)\{«„(7r)} 



7r'GSn_l 



Here the first equality follows by (|2.4|) . and the second equality follows from the bijection 

IcGC(7r)\{«„{7r)} 



5'^^'* 9 TT i-> vr' = Un{'^)Y[ceCM\fu (n)}^ ^ 5'n-i- Hencc ()2.3p follows, and we end the proof 



of Proposition [3l 

3 Proof of Theorem [1] 

Let mi, . . . , m„ and x he d x d matrices. Given a cycle c = (cr — )• ci — > C2 — )• • • • ^ c^) on [n] 
we define 

Rc{x; nil, ■ ■ ■ , fnn) = tr {xmcxXmc2 ■ ■ ■ xnic^) . 
More generally, for a permutation vr G S^, we define 

Rt,{x; nil, ... , rnn) = mi, . . .,mn). 

ceC(n) 
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For example, if n = 6 and 7r = (l^>5— 7-3^ 1)(2 6 ^> 2)(4 — )• 4), then 

Rtt{x; mi,m2,ms,m/i, ms, mg) = tr {xmixm5xms)tr {xm2xmQ)tr (xm4). 

The fohowing proposition, given in |GLM2| . is our starting point for the proof of Theorem[TJ 
Let d, /3, a be as in Introduction. 

Proposition 4. Let W Wd{(3 , o-,M.) and let si, . . . ,Sn G Sym((i). Then 

E[tr {Wsi)tT {Ws2) • • ■ tr (Wsn)] = • • • > ^n). 

Proof. See Proposition 1 in [GLM2j . See also Theorem 1 in jLMl) . □ 

Theorem [1] is a consequence of Proposition [3] and Proposition [3l For 1 < a,b < d, denote 
by Eab = the matrix unit of size d, whose (i,j)-entry is {Eab)ij = ^ai^bj- We apply 

Proposition m with Sj = {Ek^^_^k2j + /2 (1 < j < n). Since W is symmetric, we have 

tniWsj) = (Wk2j_ik2j + Wfc2jfc2j-i)/2 = Wk2j_ik2j, and therefore it follows from Proposition H 
that 

nWk,k2Wk,k,---Wk2^^,k2j 
= 2~" ^ P"^'''^ R^{a; Ek^k2 + -E'fefei; • • • > -E'fcan-lfean + -^fc2„fc2„-l)- 

From Proposition [3l in order to prove Theorem [H it is sufficient to show 

(3.1) Rnicr; Ek^ki + Ek2ki,- ■ ■ , ^fc2„-lfc2n + -^fc2nfc2n-l) = {{^ kpkq) p,qe[2n]) 

for any permutation n £ Sn- 

To show ()3.ip . let j4 = (^p5)p,gG[2n] be a symmetric matrix and let c = (c^ — > ci ^ C2 — > 
•••—)■ Cr) be a cycle. Equation (j3.ip follows from 

(3.2) tr (^(^2ci-l,2ci + ^2ci,2ci-l) • • • A{E2c,-l,2cr + ^2c,-l,2cJ) = Pc{A), 

with ^ = (o'fcpfc,)p,gg[2n]- Here the Eab = -^Ife"^ are 2n x 2n matrix units. However we may show 
(13.21) as follows: 



tr {A{E2ci-l,2ci + E2ci,2ci-l) ■ ■ ■ A{E2cr-l,2Cr + E2cr-l,2Cr)) 
2n 

^i2ril (-E'2ci-l,2ci + -E^2ci,2ci-l)jlJ2^j2i3 ■ • ■ ^i2r-2i2r-l (-E^2Cr - 1 ,2Cr + -E^2Cr ,2Cr - 1 )i2r- 1 J2r 

il j'2v J2r = l 



^J2r-2J2r-l ' 

jlvj2r 



Here the last sum is over {j2k-i,j2k) £ {(2cfc — 1,2ca;), (2cfe,2cfc — 1)} (A; = 1,2, .. . ,r). Hence 
we obtain (|3.2|) and therefore (j3.1|) . It ends the proof of Theorem [TJ 
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4 Orthogonal Weingarten functions 

We review the theory of the Weingarten function for orthogonal groups; see \CM\ IMat2j for 
details. Claims in Subsections I4.1fl4.4l are also seen in [Mac! VII-2]. 



4.1 Hyperoctahedral groups and perfect matchings 

Let Hn be the subgroup in 5*2^ generated by transpositions {2k — 1^2k^2k — l) {1 < k < n) 
and by double transpositions (2z - 1 ^ 2j — 1 — > 2i — 1) • (2i 2j 2i) {I < i < j < n). The 
group Hn is called the hyperoctahedral group. Note that \Hn\ = 2"n!. 
We embed the set A^(2n) into S2n via the mapping 

A^^n N 1 2 3 4 ••• 2n \ „ 

\m{l) tn(2) m(3) m(4) • • • m(2n) J 

where (m(l), . . . , m(2n)) is the unique sequence satisfying 

m = {{m(l),m(2)}, . . . , {m(2n - l),m(2n)}} , 

m(2fc - 1) < m(2/c) (1 < A; < n), and 1 = m(l) < m(3) < • • • < m(2n - 1). 
The m G A4{2n) are representatives of the cosets gHn of Hn in S2n- 

(4.1) S2n= U mHn. 

m£M(2n) 



4.2 Coset-types 

A partition A = (Ai, A2, . . . ) is a weakly decreasing sequence of nonnegative integers such that 
|A| := Yli>i -^i is finite. If |A| = n, we call A a partition of n and write A h n. Define the length 
£(A) of A by the number of nonzero Aj. 

Given g £ 82^ we attach a graph G{g) with vertices 1,2,... , 2n and with the edge set 

{{2A; - 1,2A;} I G [n]] U {{g{2k - l),g{2k)} \ k G [n]}. 

Each connected component of G{g) has even vertices. Let 2Ai, 2A2, . . . , 2 A; be numbers of vertices 
of components. We may suppose Ai > A2 > • • • > A;. Then the sequence A = (Ai, A2, . . . , A;) is 
a partition of n. We call the A the coset-type of (7 G S2n- 

For example, the coset-type of ( 7 | g 3 2 § I 5 ) 'S's is (2, 2). 

In general, given g.,g' G S2n-, their coset-types coincide if and only if HngHn = Hng'Hn- 
Hence we have the double coset decomposition of Hn in 52^: 

(4.2) S2n = I I Hp, where Hp = {g £ S2n \ the coset-type of g is p}. 

p\-n 

Note i/(in) = Hn and \Hp\ = (2"n!)V(2^(^)zp). Here 



(4.3) 



Zp = JJr'^'-(p)m^(p)! 
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with multiplicities mr{p) = \{i>^ \ Pi = ?"}[ of r in p. 

For g € S2n, denote by K{g) the number of connected components of G{g). Equivalently, 
K{g) is the length of the coset-type of g. Under the embedding M(2n) C S2n, we may define 
G(m) and K(m) for each m G ^A{2n). They are compatible with their definitions in Subsection 



4.3 Zonal spherical functions 

For two functions /i , /2 on S2n , their convolution fi * f2 is defined by 

(/l*/2)(5)= E /l(5(5')"')/2(5') (5G52„). 

Let "Hra be the set of all complex- valued if„-biinvariant functions on S2n- 

'Hn = {f:S2n^C\ fiCg) = f{gC) = fig) {g e ^2n, C e ^n)}. 

It is known that this is a commutative algebra under convolution, with unit 1-^^ given by 



(4.4) InM 



(2"n!)-i iigeHn, 
otherwise. 



Therefore {S2n,Hn) is a Gelfand pair in the sense of [Mac^ VII. 1]. The algebra T-Ln is called the 
Hecke algebra associated with the Gelfand pair {S2n,Hn)- 

For each A h n we define the zonal spherical function uj^ by 



2"n! 

Here x^'^ is the irreducible character of S2n associated with 2A = (2Ai, 2A2, • • • ). The uj^ (A h n) 
form a basis of T-Ln and have the property 



(4.5) UJ^*UJ^' = Sx^^^^i^^ for aU A, /i h n. 



Here /^•^ is the value of at the identity of S2n , or equivalently the dimension of the irreducible 
representation of character x^'^- We denote by io^ the value of w'^ at the double coset Hp. Note 
"^(xn) = 1 for all A h n. 

4.4 Zonal polynomials 

We now need the theory of symmetric functions. Let A be the algebra of symmetric functions 
in infinitely-many variables xi,X2,--- and with coefficients in Q. Let A = (Ai,A2,...) be a 
partition of n. We denote by px the power-sum symmetric function: 

Pa = JJpa, and pk{xi, X2, . ■ ■) = Xi + X2 -\ . 

i=l 
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Let Z\ be the zonal polynomial (or zonal symmetric function): 

(4.6) Z^ = 2^n\Y,'^-'(Ph;'oj^pPp. 

p'rn 

Here Zp is the quantity defined in ()4.3p . Alternatively, for p\~ n, 

(4.7) Pp = SE/"-p^- 

^' Ahn 

Recall that A is the algebra generated by {pr \ r > 1} and that the pr are algebraically 
independent. Let z be a complex number and let i;^^ : A — > C be the algebra homomorphism 
defined by (pz{Pr) = z for all r > 1. Then we have the specializations 

(4.8) MPp) = z"^'^ and <A.(Za) = Cx{z) := \[{z + 2j-i- 1) 

(ij)6A 

where the product ]Xj^j)g;^ stands for Hil^i nj=i' which is over all boxes of the Young diagram 
of A. It follows by and (g^D that 

(4.9) C,{z) = 2«n! 2"'^^) z^'u'^/^^^ and /^^^ = |^ E /''-pC'.(z). 
4.5 Weingarten functions 

Let z be a complex number such that C\{z) / for all Ahn. We define a function Wg'^(-; z) 
in Un by 

(4.10) Wg«(5; .) = ^^^^ ^ ^"^'^^^ 

We call it the orthogonal Weingarten function (or Weingarten function for orthogonal groups) . 

The function g H- Wg'^{g; z) is constant at each double coset Hp {p\- n). We denote by (the 
same symbol) Wg^ {p;z) its value at Hp. 



Example 1. 



Wg«((l);z) =i. 



Wg^((2); z) = , , -. Wg^((l^); z) 



"z(z + 2)(z - !)■ ' z{z + 2){z-l)' 

The list of Wg°{p; z) for |p| < 6 is seen in [CM] . 
Define the function G'^(-; z) in by 

GO(5;z) = z'^(^) (<7G52n). 
The following lemma is a key in our proof of Theorem [5J 
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Lemma 5 ( |CM] ). 

G«(.;^)*WgO(-;z) = (2"n!)2l^„. 

Here is defined in ()4.4p . 

Proof. Recall that if p is the coset-type of g, then K,{g) = i{p). From the second formula in 
.91). we have 



(4.11) g^(-;-) = |^E/''c^a(.)c.\ 



so that 

(2n)! 



Ahn 

by (ITOl) and g3|). 

On the other hand, since lim^gK^ t^+oo i~^C\{t) = 1, using the second formula in (|4.9p again, 
we may see that 

^ ' Ahn ^ ^ Ahn 

which is equal to 1 if 5 € Hn, or to zero otherwise. Hence we have 

^ ' Ahn 

This finishes the proof. □ 
4.6 Weingarten calculus for orthogonal groups 

The content in this subsection will not be used in the latter sections. We here review how the 
Weingarten function Wg*^ appears in the theory of random orthogonal matrices. 

Let 0{N) be the compact Lie group of N x N real orthogonal matrices. The group 0{N) 
is equipped with the Haar probability measure dO such that d{UiOU2) = dO for fixed Ui,U2 G 
0{N) and that J^^^-, dO = 1. 

Let O = (Ojj)jjG[Af] be a Haar-distributed orthogonal matrix. Consider a general moment 

E[Oiij,Oj2j2 ■■■Oi^jJ iii,i2,- ■ ■ ,ik,ji,j2,- ■ ■ Jk e [iV]). 

From the biinvariant property for the Haar measure, we can see immediately that E[Oj^jjOj2j2 ' ' ' ^ifcifc] 
if A; is odd. 

Proposition 6 ([UMlES])- Let ii, ■ ■ ■ ,i2n, ji, - ■ ■ , j2n be indices in [N]. Assume that N > n 
and let O = (Ojj)i,jG[Ar] be a Haar-distributed orthogonal matrix. Then we have 

m,nGA1(2n) \{p,g}em / \{p,(j}en 

Here each m € A4{2n) is regarded as a permutation in S2n- 
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In particular, using Example [H we have 

IE[Oi,jiOij2 02j3 02,j4] = j\r(-j\r_^2)(A^- 1) ^^^ ^ ^)^jihh3H " ^hh^hU " ^hh^hh) 
for N >2 and ji,j2,j3,M G [A^]- 

Remark 3. Proposition [6] was first proved in [CS] with a function Wg*^, which was implicitly 
defined via the equation of Lemma EJ The explicit expression (|4.1U|) was first given in [CM\ . 
Zinn-Justin (see also |Mat2j ) gave another expression, involving Jucys-Murphy elements. 

Remark 4. If £(A) > N then Cx{N) = 0, and therefore the definition ()4.10p does not make 
sense unless unless N > n. For z = N G {l,2,...,n — 1} we extend the definition of the 
Weingarten function by 

1 f^^ 

Ahn 

e{\)<N 

Then Wg'^ {g; N) does make sense for all g G S2n, and Proposition [6] holds true without any 
condition for N; see [CMj for details. 



5 Proof of Theorem [2] 

Let d,f3,a be as in Introduction. We also use symbols defined in Section [H Our starting point 
for the proof of Theorem [2] is the following lemma. 

Lemma 7. Let W ~ Wd{l3,cr;M.) and let si, . . . , Sn G Sym(d). Put 7 = /3 — and suppose 
7 > 0. Then 

tr {a-\si)tr (a-^ss) • • • tr {a-^^) = (-1)" J] i—,r^^^E[R^{W-'; si, . . . , 
where R-iri'] • • • ) defined in Section\^ 

Proof. We can obtain the proof in the same way to [GLMH Theorem 3] . Therefore we omit it 



here. (The assumption 7 = /3 



d+l 



> implies that the real Wishart distribution Wd^p^a has 



the density f{w;d,l3,a) given by (jl.ip . and that f{w;d,/3,a) vanishes on the boundary of $7. 
Therefore, we can apply Stokes' formula for /; see pages 298-299 in jGLMlj .) □ 

Lemma 8. Let W and 7 be as in Lemma^ Given indices ki, k2, ■ ■ ■ , /c2n from {1, . . . , d}, we 

have 



(5.1) 



' a 



k2n-lk2n 



(-1)"2 



meA^(2?i) 



n 



{p,q}&m 
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Proof. By using Lemma [71 one can prove it in the same way to the proof of Theorem [TJ Indeed, 
applying Lemma [7] with sj = {Ek^^_^^k2j + ^kijMj-i)/'^ (1 < j < "-), and using ([STID and 
Proposition [3l we see that 



= (-l)"2-" {-ir'^MRAW-'-.Ek,k2 + Ek2k,,. . . , Ek2„_,k2„ + Ek2„k2„_,)] 

ttGS„ 

= (-l)"2-"E [hf_27(W^'^'' W[2n]" • 



□ 



Suppose 7 > n — 1. Then Wg'^{g; —2j) {g € S2n) can be defined (see Subsection I4.5p . Set 

(_-,\n2n 

{2ny} ' Z^Ca(-27)- 



(5.2) Wg(5;7) = (-l)"2"WgO(5;-27) = -— (-ir2"V ^; ^ y {g) {g G S2n) 



Xhn 



We finally prove Theorem [2l Recall that the functions g i— )• K{g) and g i-^ Wg{g; z) are 
//„-biinvariant. We can rewrite (15. ID in the form 



by the coset decomposition (|4.ip . Therefore, the right hand side of (|1.4p is equal to 
(-l)"2"(2"n!)"^ ^ Wg'^(5';-27)f7V(i)V(2) ...o-V(2n-i)V(2n) 

=(2"n!)-2 (-27)'^(3)Wg«(5';-27)E 

=(2"n!)-2 (-27)'^(^?)WgO(/5"^-27)E 



VF''9'9(1)V9(2) . . . ]^V9(2n-l) Vs(2n) 



|^fcg//(l)A:g„(2) . . . -I^V(2n-l) V'(2n) 



by letting g" = g'g. Since Lemma [5] implies 



g&S2n 



2-n! iig"eHn, 
otherwise. 



the last equation equals 
(2"n!)-i Y ^ 



l^V'(l)S"(2) . . . ^'«9"(2n-l)'«9"(2n) 



Hence we have proved Theorem O 
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Remark 5. Theorem [2] holds true for any positive real number 7 such that C\{—2^) 7^ for 



Remark 6. The complex- Wishart version of Theorem [2] is obtained by Graczyk et al. [GLMlj . 
They employ a class function on Sn defined by 



where g > n — 1 is a parameter in [GLMlj . corresponding to our 7. The function Wg^(7r;7V) 
coincides with the Weingarten fucntion for the unitary group U{N), studied in [C] (see also 



6 Applications 

In this section, we give applications of Theorem [1] and Theorem [2j 
6.1 Mixed moments of traces 

Recall the symbol Rn{x; mi, . . . , m„) defined in Section [3l where x is a dx d symmetric matrix, 
nil, ■ ■ ■ , mn are d x d complex matrices, and n G Sn- For example. 



i?(i_^4_^5^i)(2^7^2)(6^6) {x] "T-i, ^2, • • • , mj) =tr (xmixm4xm5)tr (xm2xm7)tr (xnie). 

Thus Rtt{x; mi, . . . , mn) is a product of traces of the form tr {xm-i^xmi^ ■ ■ ■ xmi^). Our purpose 
in this section is to compute moments of the forms 



all A h n. 




[MN]). 



R{i^z^2^4:~^i){x; mi, m2, m-a, 7714) =tr {xmixm^xm2xm/i) 



¥.[R^{W;mi,...,mn)] 



and 



W.[R^{W-^-mi,...,mn)] 



where W ~ Wd{P,(7;M) as usual. 

First we observe a simple example. 



Example 2. We compute E[tr (H^miVFm2)]. Expanding the trace, we have 



W.[iT{WmiWm2)\= ^ {mi)k2k3{m2)kik^'^[Wkik2Wkzki]- 



ki,k2,k3,k4 



From Theorem [T] or (jl.3p , it is equal to 




:/3^tr {amiam2) H tr {am\am2) H tr (ami)tr {am2). 
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where m* is the transpose of m. In other words, 

iE[i?(i^2^i)(^; mi,m2)] = /3^i?(i^2^i) (o"; "^i, "i2)+f -R(i^2^i)(o-; "li, "i2)+^^(i^i)(2^2) (o-; mi, m2) 



This example indicates that we should deal with not only mi, . . . , m„ but also with their trans- 
poses m\, . . . , m^. 

Given a matrix m = {rriij) and a signature e G {— 1,+1}, we put 



m 



m if e = +1, 
m* if e = — 1. 



Let mi, . . . ,m„ be d x d complex matrices and let x = (xij) he a d x d real symmetric 
matrix. Given a permutation g € S2n, we define Tg{x] mi, . . . , mn) by 

d 

Tg{x;mi, . . . ,mn) = ^ (™'l)jlj2(™'2)j3j4 • • • ("'-n)j2n-l,j2n^i9(l),jg{2)^J9(3).i9(4) ' ' ' ^i9{2n-l) >i9(2n) " 

ilvj'2n = l 

In our situation, the symbol Tg is more useful than R-,^. 

Given tt € 5^, we denote by vf the permutation in 5*271 given by 7r(2j — 1) = 27r(j) — 1 and 
7f(2j) = 2j for j = 1, 2, . . . , n. Denote by the transposition (2z — 1 ^ 2i ^ 2z — 1). 

Lemma 9. For vr E and ei, . . . , e„ € {±1} we have 

Rn{x;ml\ . . . ,m^^) = Tg{x;mi, . . . ,mn) with g = I Ci • 

\j:ei = -l / 

Proof. First we will show 

(6.1) i?^(x;mi, . . . ,m„) = Ti{x;mi, . . . ,m„). 

Take a cycle c = (ci ^ C2 ^ • • • ^ — > ci) in vr. Then we see that 

r 

^ n ("^^fe ) Ja^fe - 1 Ja^fe ^is(2cfc - 1) J#{2cfe ) 

j2ci-l,i2civ:i2cr-l j'2cr k = l 
T 

~ 5^ n(™'Cfe)j2cj.-l,i2cfe^j27r{cj.)-lj2cj. 

j2ci-l J2civ:i2cr-1 :i2cr ^ = 1 

= (™Ci)j2ci-l J2ci^i2ci J2C2-1 ('^C2)j2c2-lJ2c2^i2c2 J2C3-1 ' ' ' ("^Cr )i2cr-l J2cr ^J2cr J2ci -1 

j2ci - 1 J2ci , • • • j'2cr - 1 j'2cr 

=tr {mciXmc2X ■ ■ ■ mc,.x) = Rc{x; mi, . . . , mn)- 



We obtain (|6.ip by taking the product over all cycles in tt. 
Next we will show 

(6.2) i?^(x;mi, ...,m\,.. . ,mn) = TQi{x;mi, . . . ,m„). 
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We have 

n 

jl,---,j2n k = l 



Letting = jQ(k) for all /c = 1, 2, . . . , 2n, it is equal to 



5^ IT ^^'^•^■?Ci(2fc-l) '^'Ci(2k) '^■'4{2fc-l) '•'ir(2fc) 

^ "-'2i-l'-'2i J#(2i-1) 'J7r(2i) ±1^ '^'■J2k-l'J2k ■Jw(2k-l)'-'n{2k) 

=T^{x;mi,. . . ,m-, . . . ,mn). 



Therefore, ()6.2p follows by ()6.ip . Now the result can be obtained from ()6.ip and ()6.2p . □ 
Example 3. Consider 

tr {xmixm\xm\xm2)iT (xmaxm^jtr (xme), 

which is equal to Rnix; m'^ , . . . , m^^) with 

7r=(1^4^5^2^ 1)(3 ^ 7 ^ 3)(6 ^ 6) e 5/, 
(ei,...,e7) =(+!,+!, +1,-1,-1, +1,-1). 

It coincides with Tg{x;mi, . . . ,mj), where g = CiCbCr^i i-e-j 

5 = (7 ^ 8 ^ 7)(9 ^ 10 ^ 9)(13 ^ 14 ^ 13)(1 ^ 7 ^ 9 ^ 3 ^ 1)(5 ^ 13 ^ 5)(11 ^ 11). 

Lemma 10. The function S2n ^ g ^ Tg{x] mi, . . . , m„) is right Hn-invariant: 

Tgc_{x]mi, . . . ,mn) =Tg(x;mi,...,m„) for all Q E i7„ and g G S2n- 

Proof. It is enough to check for C = (2i - 1 ^ 2i ^ 2i — 1) and (2i — 1 ^ 2j — 1 ^ 2f - l)(2i ^• 
2j — > 2i) because is generated by them. However, it is clear. □ 

The moment of the form '¥\R.,^{W^^]nv^ , . . . , rrf^)] may be given by W(Tg{W^^]mi, . . . , m„)] 
with some g £ S2n- Hence we now compute the moments 'K\Tg{W^^]mi, . . . , m„)]. First of all, 
we note that the formulas in Theorem [T] and Theorem [2] can be expressed in the forms 

(6.3) nWk.k, ■ ■ ■ Wk,„,,,,J =2-"(2"n!)"i Yl (2/3)''^'^^fc,(,),fc,(2) ' ' ' ^/^,(2„-i),/^,(2„) > 

gG52n 

(6.4) E[VF''i''2 • • • "H/^2n-l,A:2n] ^(-2«„!)-l ^ Wg(^;-)/)o-''9(l)''=9(2) . . .o-'=9(2n-l)'*:9(2n)_ 

9e'S'2„ 
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Theorem 3. Let W ~ Wd(/3, fT;M) and let 7 be as in Theorem Let rn-i,...,m„ be d x d 
matrices and let g € S2n- Then 

nTg{W;mi,...,mn)]=2-^ ^ (2/3r(^"'")r„(a; mi, . . . , m^, 

ne>1(2n) 

E[rg(W^-i;mi,...,m„)] = ^ Wg((7-'n; 7)r„((T-i; mi, . . . , m„). 

n6A^(2n) 

Proof. Using (|6.3p (or Theorem [T]), 
E[r3(Ty;mi,...,m„)] 



ilvj'2n \fc=l / 



2n)J 



{2n) 



ilvj2n \fc=l / 3'eS2,i 

and, letting /i = (7(7', 

n 

=2-"(2"n!)-i {2(3)^^-'^^ Il(^k)n.-.,n.^M^.-.yM2., 

h£S2n jl,---,j2n k=l 

=2-"(2"n!)-i m^^'''''^n{a;mi,...,mn) 
=2-" (2/3r^^"'"^r„(a;mi,...,mO. 

neX(2n) 

Here the last equality follows from Lemma fTOl and ()4.ip . Thus the first formula has been proved. 
The same applies to the second formula. □ 

It follows from Lemma [9] and Theorem [3] that, for vr E 5'„ and (ei, . . . ,e„) G { — 1,+!}"', 

(6.5) nRAW;m\\...,m'-)]=2-^ ^ (2/3r(f "")Tn(a; mi, . . . , m„), 

neM(2n) 

(6.6) ^[R^{W-m\\...,m'-)]= Y Wg((7-in; 7)r„(a-^ mi, . . . , m„), 

nGA1(2n) 

where 3 is as in Lemma [9j We remark that (|6.5p is equivalent to |GLM2[ Corollary 14]. 
6.2 Averages of invariant polynomials 

Given a partition A of n, we define two functions Zx and px on = Syin^{d) by 

Zx{x) = Zx{ai,a2,.. .,0^,0,0,...) and px{x) =px{ai,a2, . . . , a^, 0, 0, . . . ), 
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where oi, . . . , a^^ are eigenvalues of a; G fi, and Zx,px are symmetric functions defined in Sub- 
section In particular, we have 

i=l r>l 

where mj.(A) is the multiplicity of r in A. From (j4.6p and ()4.7p . we have 

(6.7) Z, = 2"n!5]2-^(^)^;ia.>, and = ^ f^u^^^Z^. 

p\-n ^ Xrn 

Recall Cx{z) = J^^^ + 2j ~ i — 1). The following theorem, derived from Theorem [1] and 

Theorem [2l is exactly the real case of Proposition 5 and 6 in [LMlj . 

Theorem 4. Let W ~ VFd(/3,(T;M) and let 7 he as in Theorem\^ For a partition A of n, 

nZx{W)] =2-"CA(2/3)ZA(a). 
nZx{W-^)] =(-l)"2"CA(-27)-iZA(a-i). 

Proof. First of all, we note that 

Pp{x) = Tg{x; Id, . . . , Id) 
«, ' 

n 

for a permutation g in S2n of coset-type p and for a matrix a; in fi. Indeed, since the function 
S2n ^ 9 '-^ Tg{x; Id, ■ ■ ■ , Id) is i/„-biinvariant, the image depends only on the coset-type. If vr is 
a permutation in Sn of cycle-type p, then n is of coset-type p, and therefore rg(x; Id, ■ ■ ■ , Id) = 
(x; Id, ... , Id) = R^ix; Id, ■ ■ ■ , Id) = Pp{x) by Lemma [3 

From the first formula in (j6.7p and the double decomposition (|4.2p . we have 
E[Za(T4^)] =2-n! J] 2"^('')z;i.;^E[p,(IF)] 

phn 

=2"n! j;2-^('')z;i-^ ^ ..^(g)E[T,(H^;/d,...,/d)] 
= (2"n!)-i ^ a;^(5)E[r,(I^;/d,...,/d)]. 

It follows from Theorem [3] that 

E[Zx{W)]={2^n\)-^ u\g)2-^ ^ (2/3)«(f "V)T,,(a; J^, . . . , 1^) 

=(2"n!)-22-" ((.;^*G«(.;2/3))(5'))T,,(^;/d,...,/d). 
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Since u we have 

E[ZA(t^)] = (2"n!)-i2-CA(2/3) J] L.^(5')Tg. (a; /,,..., 



Since 

/on |-\2 

by the first formula in (|6.7|) . our first result follows. The proof of our second result is similar. □ 

The following is equivalent to the real case of |LMH Theorem 2]. 
Corollary 5. Let W ~ Wd{(3,cr;M.) and let 7 be as in Theorem\^ For a partition /i of n, 

^ '' p^n \ Ahn / 

^ p\-n \ Ahn / 

Proof. They follow from Theorem S] and ()6.7p . □ 
Corollary 6. Let W ~ Wd{/3,a;M.) and let 7 be as in Theorem\^ Then 

E[(trT^)"]=E-/5^(^)pp(a), 

p^n ^ 

E[(tr W~^T\ = 2"-^(^) -Wgip; i)pp{a-^). 

phn ' 

Proof. The first result follows by letting fi = (1") in Corollary [S] and by using the second formula 
in ([321) • The second one also fohows by (j4.10p . □ 

6.3 Examples for low degrees 

If we apply our theorems (Theorem [1]- Corollary [6]) , we can compute various Wishart moments 
for low degrees n = 1,2,3,4 easily; see Appendix. 
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A Appendix: Examples for low degrees 

We give explicit examples of our theorems obtained in the present papers. Let W ~ Wd{(3, a; ] 
and set 7 = /3 — as usual. Let mi,m2, • • • he d x d matrices. 

Degree 1 

Suppose 7 > 0. It follows from Theorem [1] and Theorem [2] that 
(A.l) E[Wij] = (3aij and E[W'^] = ^a'^ 



for i < i, j < d. It is immediate to see that 
(A.2) K[W] =ficj, 



(A.3) 



E[tr(H^mi)] =/3tr(cTmi), 



E[tr(W^"imi)] =7"Hr [a'^mi) 



Degree 2 

Suppose 7 > but 7^1 (see Remark [S]) . Prom (|5.2p and Example [TJ 



Wg({{l,2},{3,4}};7) 
Wi({{l, 3}, {2, 4}}; 7) = Wg({{l, 4}, {2, 3}}; 7) 
It follows from Theorem [1] and Theorem [2] that 



27 - 1 



"7(7-l)(27 + l)' 
1 

>(7-l)(27 + l)' 



(A. 4) E[pyfc^jt2^fc3fc4] 0'fciA:2^fe3fc4 + "^(<7fclfc3'7fc2fc4 + (^fci fc4 0'A:2 fca ) ) 



(A.5) E[iy''i''2|^''3fc4 



7(7-l)(27 + l) 



^2^1 _ \'^(jklk2^kzk4, _j_ ^kik-i ^k2k4, _j_ ^kik4^k2k3 



for {ki,k2,kz,ki) G [d]^. 

The average for the (i, j)-entry of VF^ is 



E 



Lfc=i 



k=l 



k=l 



and the average for the (i,j)-entry of ^ is 

d 



E 



.k=l 



1 



7(7-l)(27 + l) 
1 



(27 - 1) ^a^V^J' + ^((T^V'^-'' + a'^a'''') 



k=l 



k=l 



7(7-l)(27 + l) 



-(27(0.,+tr(a-iK^-) 
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Therefore 

(A.6) ElW-'] = (0' + ^y' + ^{tva)a, 

As we saw in Example [21 

(A. 8) E[tr (WmiWm2)] = /3^tr {amiam2) + — tr {am\(Tm2) + — tr ((Tmi)tr (am2), 



and in a similar way we have 



1 



(A.9) E[tr {W-'miW-'m2)] =—, , (27 - l)tr {<j- ' mia- ' 7112) 

7(7- 1)(27 + 1) L 

+ tr {(7~'^m\a^^m2) + tr {a~^mi)ti {a~^m2) 



Moreover 



(A. 10) E[tr (l^mi)tr (14/^m2)] =/3^tr ((Tmi)tr ((77712) + ^tr (177771(77772) + ^tr ((7777*cJ7772), 

(A.ll) E[tr (l^-^777i)tr (1^-^7772)] =— ^7 [(27 - l)tr (cr-i777i)tr (ct" V2) 

7(7 - 1)(27 + 1) L 



+ tr ((7 mia 7772) + tr ((7 niiG 7772) 



Degree 3 



Suppose 7 > but 7 7^ 1, 2. From ([O]) and a hst in fCM] (see also [CS]), the Wg(/9; 7) (p h 3) 
are given by 

where 

^3(7) = 7(7 - 1)(7 - 2)(7 + 1)(27 + !)• 
It follows from Theorem [1] and Theorem [2] that 

(A.12) E[PFfe,fc,PFfc3fc,Tyfc,fcJ 

=/3 0-kik2^k3k4,0-k^kG + -^(c^fcifc3<7fc2A:4<7A:5fc6 + '^kik4,0'k2k30'ksk6 + O-kiks'^kike'^kski 
+ cr kikfiCr k2k5(^k-iki + Cr kiki^k-iks^^kika + kik2^kak(i'^ kiks) 

+ ^(Crfcifc4'^fc2fc5'^fc3fc6 + ^kxka<^k2k5(^ kike + (^kik4(^k2ke'^k3k5 + <^kihi(^ k2ke<7kAkr, 
+ cr kikfi<y k2k3Cr kik^ + kik^cr k2k3(^kAkfi + <y kikfi<^k2kA(^k'ik5 + <y kik5(^k2kA<^k-ika) 
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and 

(A.13) E[W''^''^W''^'''^W''^'"^ 



=^3(7)"^ (27^ - 37 - l)f7''-'-0-— 0- 

_)_ (-y _ I'^^fjkik-i^k'zki^kskf} _|_ ^kik^ ^k2kz ^k^kg, _j_ ^kikz, ^k2kfi ^k^k^ 

_j_ ^kik^ ^k2kz, ^k-iki _|_ ^kik2 ^ks,k^ ^k^k^ _|_ ^kik2 ^k^k^ ^kik^^ 
_j_ (^^kiki ^k2k^ ^k^kfi _|_ ^kiks ^k2kr, ^k^ka _|_ ^kiki ^k2k(i ^k^k^ _j_ ^kikz ^k2k(i ^k^k^ 

_|_ ^kikfi ^k2k3 ^kiks _|_ ^kiks ^k2k3 ^k^kg _j_ ^kika ^k2k4 ^k-^ks _|_ ^kiks ^k2k4, ^kikg') 

From Corollary [5] we have 
(A.14) 

=y (3))P(3) (a) + ^^(//, (2, l))p(2,i) (a) + ^^(/x, (l3))p(i3) (a)) , 

(A.15) 

=f (3))p(3)(a-i) + (2, 1))P(2,1)(CT-1) + (I3))p(i3)(c7-i)) 

for each ^ H 3, where 

Mf^,P) = J^CA(2/3)/2^a.;i^^^ and = -8 J] Ca(-27)-V 



-1 ^2A A A 



Ah3 AI-3 



We compute the matrices A = {A{fi, p))^ p\-^ and B = p))^^pi-3. Here indices of rows and 

columns of the matrices are labeled by (3), (2, 1), (1'^) in order. By using results in [Maci VII. 2], 
we have 

1 11^ 



Z ■= (w^)a,mI-3 = I -4 6 1 

-I h 

Since /^^ coincides with the number of standard Young tableaux of shape 2 A (see, e.g., [Sajl ), 
we may have 

/2(3) = /(6) ^ 1^ J2(2,l) ^ ^(4,2) ^ ^2{l3) ^ ^(23) ^ 

From the definition of Cx{z), it is immediate to see 

C7(3)(z) = z(z + 2)(z + 4), q2,i)(z) = z(z + 2)(z-l), and C^^.){z) = z{z - l){z - 2). 

Now, letting F := diag(/2(3), /2(2,i) j2(i3)) diag(C(3)(z), C(2,i)(^), C(i3)(z)), we can 

calculate 

/i|/3(2/32 + 3^ + 2) f/3(2/3 + l) f;9' 
^= -Z*-F-(:7(2/3) •Z= I f/3(2/3 + l) f /3(2/32 + /3 + 2) f/S^ 

^(3 ^/32 15/33 



^[Sa] B. E. Sagan, The symmetric group. Representations, combinatorial algorithms, and symmetric functions, 
second ed.. Graduate Texts in Mathematics, 203. Springer- Verlag, New York, 2001. 
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and 



''il' f7 15 



S = -8Z* • F • C(-27)-i • Z = ( f 7 5(72-7 + 1) 15(7-1) 



^^3(7) 



15 15(7 - 1) 15(27^ -37-1) 



Hence 



(A.16) E[p(3)(W^)] =^/3(2/32 + 3/3 + 2)p^,){a) + ^/3(2/3 + l)p(2,i)(^) + ^/3p{i3)(a), 

(A.17) E[p(2,i)(H^)] =/3(2/3 + l)p(3)(^) + ^/3(2/32 + /3 + 2)p(2,i)(a) + ^/32p(i3)(^7), 
(A.18) ]E[p(i3)(VF)] =2/3p(3)(a) + 3/32p(2,i)(a) + /3='p(i3)(a), 

and 

(A.19) E[p(3)(T^ )]- ,(^_i)(,_2)(7 + l)(27 + l) ' 

f A 2m Fr« fM/^i^l _ 47P(3)(^-^) + 2(7^ - 7 + l)P(2,i)(^-^) + (7 - l)P(i3)(^-^) 
(A.20) E[p(2,)(H^ )]- 7(7-l)(7-2)(7 + l)(27 + l) ' 

rA2n Ffr, (w-^^^ - ^P(3)('^"') + ^(7 - i)p(2,i)(^-^) + (27^ - 37 - i)P(i3)(a-^) 

(A.21) E[p(,3)(M^ )]- 7(7-l)(7-2)(7 + l)(27 + l) " 

We remark that those formulas for E[p^(VF)] (/i h 3) are seen in [LMTI equation (37)]. 

Degree 4 and higher degrees 

First we note that, when n = 4, the sums in Theorem [H [2] and [3] are over |7W(8)| = 7-5-3-1 = 105 
terms. 

Consider Corollary [5] for any degree n. As we did in the degree 3 case, we can apply it 
to any degree n. The f"^^ may be computed by the well-known hook formula, see e.g. [Sa, 
Theorem 3.10.2], and the C\{z) may be done easily by the definition ()4.8p . The a;'*' are the most 
complicated among quantities appearing in Corollary [5] but we can know their explicit values 
from the table of zonal polynomials in 

In closing, we give the explicit expressions of Corollary [6] for n = 4. Its first formula is given 

(A.22) E[(trW)^] = 6/3p(4)(cT) + 8/32p(3^i)(a) + 3/32p(22)(a) + 6/33p(2^i2)(a) + /34p(i4)(a). 
Suppose 7 > but 7 ^ |, 1, 2, 3. Put 

n4(7) = 7(7 - 1)(7 - 2)(7 - 3)(27 - 1)(7 + 1)(27 + 1)(27 + 3), 



[PJ] A. M. Parkhurst and A. T. James, Zonal polynomials of order 1 through 12, Selected Tables in Mathe- 
matical Statistics (1974), vol. 2, 199-388. 
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which is non-zero. From ([521) and a Hst in [CM] ( see also |CS| ). we have the exphcit values 
Wg((4);7) =^^, Wg((3,l);7) -^^^^"^^ 



^4(7) ' ' ' ^^4(7) 

Wg((2=);7)=^2^^. Wg((2.1^);-,)=i2!^i?2!±^, 

^4(7) «4(7) 

«4(7) 

Hence the second formula of Corollary [6] at n = 4 is given 

(A.23) n4(7) • E[(trTy-i)4] =48(57 " 3)P(4)(ct"') + 1287(7 " 2)P(3,i) (^"') 

4-12(272 -57 + 9)p(22)(a-i) 

+ 12(473 - 1272 + 37 + 3)p(2,i2)(a-i) 

+ (7 + 1)(27 - 3)(472 - 127 + l)p(i4)(a-i). 



